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Abstract. Event-by-event fluctuations and correlations in azimuthal angle are currently 
widely investigated in various experiments. In this paper the $ measure (earlier used in 
experiments to evaluate fluctuations in transverse momentum) is now applied to azimuthal 
angle <f>. Properties of this function are investigated through fast generators and with 
complex models such as Pythia, Hijing, and UrQMD. Preliminary results of NA49 on $0 are 
also presented. The system size dependence (p+p, C+C, Si+Si and 6 centralities of Pb+Pb) 
at the highest SPS energy (158^4 GeV) is shown, as well as the energy dependence (20^4 - 158A 
GeV) for the 7.2% most central Pb+Pb interactions. 
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1. $0 measure 

1.1. Motivation 

The reasons to perform event-by-event azimuthal angle fluctuations were to investigate plasma 
instabilities [1], search for the Critical Point and Onset of Deconfinement, and measure 
fluctuations of elliptic flow [21 [3] . In this work the measure [3] is chosen as it is a strongly 
intensive measure of fluctuations (does not depend on volume and on volume fluctuations). The 
$ measure was already successfully used by NA49 to study average pt fluctuations |6] and 
charge fluctuations [7|. There are several other effects that influence the values of this measure, 
such as resonance decays, momentum conservation, flow, and quantum statistics. The <3? measure 
is not well suited to the analysis of these effects which represent a background for the present 
study. 

1.2. $0 definition 

Let <j) be a particle's azimuthal production angle. One can define a single-particle variable 
za = <p — 4>, where (j) is an average over the single-particle inclusive azimuthal angle 
distribution. Let us also define an event variable = Y^i=i{4>i ~ 4>)i where the 

summation runs over particles in a given event. Then the <E>0 measure is defined as: 

where (...) represents averaging over events. When particles are emitted by a number (N s ) of 
identical sources, which are independent of each other and P(N S ) is the distribution of this 
number, then ^(iV,) is independent of (N s ) (an intensive measure) and its distribution P(N S ) 
(a strongly intensive measure). Additionally, if particles are produced independently (no inter- 
particle correlations), <3?0 equals zero. 

2. Model studies 

General properties of the <3?^ measure were already studied using analytical calculations [8]. 
Fluctuations due to elliptic flow, quantum statistics and resonance decays in the hadron gas 
were evaluated. In this paper we present model studies that were conducted to investigate in 
detail the behavior of $</,. Several toy models were employed to test separately the influence of 
various physics effects. Full event generators such as UrQMD [9j, Pythia [lUj . and Hijing [TT] 
were used as a direct reference to experimental results- 



see conference slides for Pythia and Hijing results 



2.1. Toy models of elliptic flow and momentum conservation 

To simulate elliptic flow events were generated with azimuthal angle distribution of particles 
following p((p) = 1+2^2 cos (2(0— <t>it))- For each event the reaction plane angle <j)R was generated 
from a flat distribution in azimuthal angle, and the multiplicity N of particles in an event was 
taken from a Negative Binomial distribution with given (iV) and dispersion Dn ~ 0.5- {N). The 
value of V2 was a simulation parameter, constant for each simulation series (Fig. [1]) or varying 
from event to event according to a Gaussian distribution with a V2 (Fig. [2]). 
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Figure 1. $0 as a function of t>2 for 
elliptic flow simulation with constant v 2 . Lines 
correspond to analytical formulas from [8], 
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Figure 2. Difference between $^ values for 
constant and fluctuating v-i as a function of 
relative size of i>2 fluctuations. 



The elliptic flow effect results in positive values which are increasing with increasing V2- 
Moreover, event-by-event fluctuations in t>2 cause an additional increase of <&0 values. 

To simulate global momentum conservation (Fig. [3D, each particle had its transverse 

—PT 

momentum px drawn from the distribution P{pt) ~ Pt£^~ , where T = 200 MeV, and its 
azimuthal angle from a flat distribution. Afterwards, for every particle p x and p y was modified: 

EN r-^JV 
Px — Px Xl ; and p y = p y — f$ Vt to obey momentum conservation in the whole event 

(iV is the event multiplicity generated from a Negative Binomial distribution). Results show 
negative values (anti-correlation) which weakly depend on multiplicity. 

2.2. UrQMD - a full event generator 

A simulation using UrQMD 3.3 was performed for p+p collisions in a wide energy range (SPS, 
RHIC, LHC) (Fig. |4|). Default parameters of the generator were used, which for this version 
means simulation of hard processes (for energies -^/s > 10 GeV), but no hydrodynamics stage. 

For all energies and charge combinations ^ values are negative (domination of anti- 
correlations). Their magnitude increases for negatives and decreases for positives and all charged. 
At LHC energies the values become charge independent. 

3. NA49 experiment 

The NA49 experiment [32] at the CERN SPS operated between 1994 and 2002. The registered 
data include p+p, C+C, Si+Si, and Pb+Pb interactions at a center of mass energy of 6.3 - 
17.3 GeV per N+N pair. It is a fixed target hadron spectrometer (four TPCs, two of them, 
VTPC-1 and VTPC-2, inside a magnetic field) with precise determination of the collision 
centrality (Forward Calorimeter measuring energy of projectile spectators). 
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Figure 3. $0 as a function of mean 
multiplicity for simulation with momentum 
conservation. 
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Figure 4. $^ as a function of energy for 
p+p collisions generated using the UrQMD 3.3 
model. 



3.1. Detector 

The NA49 detector has a limited azimuthal angle acceptance; the acceptance losses are 
concentrated mainly around the up-down regions. Due to the magnetic field in VTPC-1/2 
particles of opposite charge are almost completely separated in the MTPCs (Fig. [5]); the MTPCs 
cover mainly the forward rapidity region. As the detector is left-right symmetric, acceptance 
for positively and negatively charged particles is the same, provided the azimuthal angle for 
one charge is reflected. To allow a quantitative comparison of ^ for positively and negatively 
charged particles we rotate particles of one charge (Fig. [6]). 




Figure 5. Tracks registered in NA49; four 
TPCs with marked regions where particles of 
given charge are dominant. 




Figure 6. Examples of <f> distributions at 
forward rapidity for positively and negatively 
charged particles before (upper plots) and after 
(lower plots) the shift. 



3.2. Preliminary results 

Figure [7] shows the energy dependence of for the 7.2% most central Pb+Pb interactions. 
There is no significant energy dependence, however negatively charged particles show positive 
<3?^ values. The results from UrQMD, obtained with the same kinematic restrictions, show values 



consistent with zero. In the simulation the limited NA49 acceptance causes a major reduction 
in the absolute values of the measure (Figs. [9land fTUl in Fig. 21 showing UrQMD results for 
p+p data, the complete rapidity and azimuthal angle range was used). The main effect is due 
to the rapidity cut. The additional azimuthal angle restriction results in a <3?0 value consistent 
with zero. 

Figure[8]presents the system size dependence of $0 at 158A GeV (top SPS energy). Significant 
positive values of $0 are observed with a maximum for peripheral Pb+Pb interactions. This 
result is qualitatively similar to the results of multiplicity [T3] and average px IS] fluctuations in 
NA49. Further studies are planned to understand these intriguing effects. 
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Figure 7. $^ as a function of energy. 
NA49 data compared to UrQMD 1.3 with 
the same kinematic cuts. See [B] for 
details of kinematic cuts and precise azimuthal 
acceptance description. 
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Figure 8. $0 as a function of system size. 
See [5] for details of kinematic cuts and precise 
azimuthal acceptance description. 
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Figure 9. $^ (for pos. charged particles) 
as a function of energy for Pb+Pb collisions 
generated using UrQMD 1.3 with different 
kinematic restrictions. 
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Figure 10. (for neg. charged particles) 
as a function of energy for Pb+Pb collisions 
generated using UrQMD 1.3 with different 
kinematic restrictions. 



4. Summary 

Event-by-event azimuthal fluctuations were analyzed using the <J>^ measure, for models and for 
NA49 data. The properties of ^ were studied using toy models and UrQMD. Elliptic flow 
results in positive values of which increase with the magnitude of flow and are sensitive 
to its fluctuations. Contrary to flow, momentum conservation produces negative values. 
The UrQMD model shows negative values with a weak energy dependence. The NA49 
measurements for central Pb+Pb collisions showed weak energy dependence and positive 3>0 
values for negatively charged particles. On the other hand, a significant system size dependence 
is observed at the highest SPS energy. This result is qualitatively similar to the behavior found 
for multiplicity [T3] and transverse momentum [5j fluctuations. 
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